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ABSTRACT 

Pitch angles p of the large-scale magnetic fields B of spiral galaxies have previously been inferred 
from observations to be systematically larger in magnitude than predicted by standard mean-field 
dynamo theory. This discrepancy is more pronounced if dynamo growth has saturated, which is 
reasonable to assume given that such fields are generally inferred to be close to energy equipartition 
with the interstellar turbulence. This “pitch angle problem” is explored using local numerical mean- 
field dynamo solutions as well as asymptotic analytical solutions. It is first shown that solutions in 
the saturated or kinematic regimes depend on only five dynamo parameters, two of which are tightly 
constrained by observations of galaxy rotation curves. The remaining 3-dimensional (dimensionless) 
parameter space can be constrained to some extent using theoretical arguments. Predicted values of 
IpI can be as large as ^ 40°, which is similar to the largest values inferred from observations, but only 
for a small and non-standard region of parameter space. We argue, based on independent evidence, 
that such non-standard parameter values are plausible. However, these values are located toward the 
boundary of the allowed parameter space, suggesting that additional physical effects may need to be 
incorporated. We therefore suggest possible directions for extending the basic model considered. 
Keywords: dynamo - galaxies: magnetic fields - galaxies: spiral - magnetic fields 


1. INTRODUCTION 

The pitch angle p of the regular (i.e. large-scale) mag¬ 
netic field may be the most important observable with 
which to test galactic dynamo theory. The quantity p 
measures how tightly_woimd is the large-scale field, and 
is given by tanp = Br/B^ where r is the galactocen- 
tric radius and (j) is the azimuthal coordinate, oriented 
along the direction of the rotational velocity. By con¬ 
vention, p is defined to be between —90° and 90°. A 
determination of p from polarized radio emission gener¬ 
ally involves some amount of inference, but p is closer to 
^direct observable than the strength of the mean-field 
B. Note that a larger magnitude of p implies more open 
magnetic field lines, p = 0 corresponds to a circular field, 
and p < 0 corresponds to a trailing magnetic spiral, p is 
generally observed to have a negative sign, in agreement 
with the predictions of mean-field dynamo theory for a 
differentially rotating galactic disk, with angular velocity 
decreasing outward. 

Observed magnetic pitch angles have magnitudes 
larger than those predicted by ki nematic mean-field dy - 
namo theory. In the data set of IVan Eck et al.l (1201511 . 
the mean value of p is —25°, and —p may be as small as 
8° or as large as 48° (excluding uncertainties), whereas 
kinematic theory generally predicts —p < 20°. More¬ 
over, this mismatch between theory and observation 
worsens in the nonlinear regime, when the fi eld ap¬ 
proaches energy equipart ition with turbulence (|Elstnerl 
I2005t IVan Eck et al.l[20T^ , and is predicted to then have 
a reduced value of —p. As local galaxies have large- 
scale fields near equipartition, it is reasonable to in¬ 
fer that they have already reached the saturated state. 
This adds up to the existence of a “pitch angle prob- 
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lem.” This problem has been addressed in the past by 
making the disc thinner, the shear smaller, assuming 
that t he kinetic a effect bec omes enhanced in spiral 
arms (iChaman dv et a lll2014bll . incorporating mean ra- 
dial fl^ s (|Moss et al.| 1200011 . or invoking spiral shocks 
(|Van Eck et al.l 1201511 . Although these proposed reme¬ 
dies may all play a role, none is likely to be sufficient 
on its own; a more natural and universal solution to the 
problem would be preferred. 

The goal of the present paper is a systematic study of 
the magnetic pitch angles predicted by standard mean- 
field galactic dynamo theory. Unlike most previous 
works, we make no a priori judgements about the dy¬ 
namo parameters; rather, we focus on mapping out the 
available parameter space and exploring the solutions 
that may lead to values of p that are in good agreement 
with observations. In Section [2] we present the dynamo 
solutions used and summarize the underlying theory. We 
then discuss important constraints on the model stem¬ 
ming from theoretical considerations in Section [3] This 
is followed by our main results in Section |4l namely the 
predictions for the magnetic pitch angle from numeri¬ 
cal as well as analytical solutions over a wide range of 
parameter values. In Section [S] we assess, citing inde¬ 
pendent evidence, the plausibility of various parameter 
values, and in Section [5] we offer several ideas for ex¬ 
tending the basic model to include additional physical 
ingredients that may have important effects on the mag¬ 
netic pitch angle. Einally, we summarize and present our 
conclusions in Section [71 

2. THEORETICAL BACKGROUND 

The evolution of the mean magnetic field is governed 
by the mean-field induction equation 

dB/dt = V X {U xB + S - T]V xB) , (1) 
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where overbars represent mean (large-scale) quanti¬ 
ties. For statistically isotropic turbulence, the first or¬ 
der smoothing closure approximation (FOSA) gives the 
mean electromotive force as 


initial values <C Beq (the solution for the mean field in 
the saturated regime is independent of the seed field so 
long as it is sufficiently small). Fo r numerical solutions, 
following iChamandv et al.l (l2014all we adopt 


S = aB-l3V xB, (2) 

where a and j3 depend on the statistical properties of the 
random (small-scale) velocity u and random magnetic 
field b. Under the same circumstances, the tur bulent 
diffusivity is given by /3 = ^ where u=\/u^ and 

is the correlation time of the random flow. Combining 
equations © and @, and neglecting t] compared with 
/3, we obtain the standard dynamo equation 

dB/dt = V X (U xB + aB-PV xB). (3) 

We adopt the dynamical quenching formalism, 
whereby the quenching of the mean-field dynamo can 
be understood to arise from the additional requirement 
of magnetic helicity balance. In this paradigm, the 
backreaction of the Lorentz force onto the dynamo is 


modeled bv including magnetic as well as kinetic terms 

in a (iPouQuet et al.1119761 IKleeorin & RuzmaikinI 

1982 

Gruzinov & Diamondl 19941 Blackman & Fieldl 

2000 

Blackman & Brandenburg! 120021 IRadler et al.l 

2003 

Brandenburg & Subramanianll2005ri . Thus, a = ak-fam. 


with 


ttk =-j'TcM • V X M, = gUMA • V X ma, (4) 

where ma = bly/Airp is the Alfven velocity associated 
with the random magnetic field and p is the gas den¬ 
sity. Because Om is closely related to the mean small- 
scale magnetic helicity density, magnetic helicity balance 
can be invoked to der ive an evolution equation for a^i 
(|Shukurov et al.ir2006[l . 

da^/dt = -(2/37^2) (5 • B/B^^) -VT, (5) 

where Bgq = u^/Airp is the equipartition field strength, 
I is the correlation scale of the random flow, and we 
have neglected an Ohmic dissipation term because it is 
negligible compared to the term involving the flux den¬ 
sity ! F of Qm. The latter is given by a sum of advec- 
tive (i Subramani an fc Brandenburg 12006 ) and diffusive 
(jBrandenburg. Candclaresi. fc Chatterieell2009l l terms 

:F = [/Om - KVOm, (6) 


with K the turbulent diffusivity of am- We note 
that other terms may be importan t in equation m 
(iBrande nburg fc SubramanianI l 2005f) or in equation 
(151) (ISubramanian fc Brandenburg) ^06t rvTshniadl 2012 t 


Vishniac fc Shapovalovll20l4lEbrahimi fc Bhattacharieel 

201411 . but they are less certain and so are left for con¬ 
sideration in a future study. 

We assume axisymmetry, adopt cylindrical (r, (j), z) 
coordinates, and make the local ‘slab’ approximation. 
That is we neglect radial as_compared to vertical c^iva- 
^es, e.g. dBz/dr dBrjdz, and assume B^ <C 
Bj.. Generally, this is justified if the local galacto- 
centric r adius greatly exceeds the disk half-thickness 
r :$> h (jRuzmaikin et al.l 119881 1. Th^resulting_ equa¬ 
tions are solved numerically to obtain Br{t,z), B^(t,z) 
Md ain_(t, z). We adopt vacuum boundary conditions 
Br = B(j, = 0 as well as B^a^dz^ at z = ±/i; we 
also set Om = 0 at t = 0 , while B^ and B^ are set to small 


Beq = Bq exp 



where Bg is the strength of the equipartition field at the 
midplane. Admittedly, exponential profiles may provide 
better fits to more observations or direct n umerical sim¬ 
ulatio ns than Gaussian profiles (see, e.g. iDumke et al.l 
119951 where both exponentials and Gaussians are used 
to model synchrotron emission in edge-on galaxies). In 
any case, the solution for B turns out to be remark¬ 
ably insensitive to the profile chosen, with vertically av¬ 
eraged p in the saturated state differing by < 1% if 
Beq = Bq exp(—|z|//i) or even Beq = Bq is us^ instead. 

The mean velocity field is taken as (0, rU, Uz), where 
n is the magnitude of the angular velocity about the 
galactic center (fl = Uz where hat denotes unit vec¬ 
tor). The observable quantities that are parameters of 
the model are the disk half-thickness h, I, u, 0, the shear 
parameter q = — 9 In 0/9 In r {q = 1 for a flat rotation 
curve), p_j_K and the amplitude Uq of the mean vertical 
velocity Uz- For the latter, we adopt the form 

Uz = Uoz/h. 

Of these, O and q are typically well-constrained by obser¬ 
vation. Further^ as_shown below, the magnetic pitch an¬ 
gle p = tan“^(Br/B 0 ) is independent of p and k. More¬ 
over, kinematic or steady state solutions can be param¬ 
eterized in terms of three ratios of the four remaining 
parameters /i, Z, u, and t/g. For our purposes, it is con¬ 
venient to consider the parameterization 

H = h/l, T = 1/it, V = Uqju. (7) 

Here B is the dimensionless half-thickness of the disk in 
correlation lengths, r is the eddy turnover time, and V 
is the dimensionless mean outflow velocity in terms of 
the rms turbulent velocity. As we shall see, r usually 
enters the expressions as the inverse Rossby number Ur 
(also equal to half the Goriolis number), which is dimen¬ 
sionless, and it is this quantity that we generally work 
with. However, it can be interesting to consider r sepa¬ 
rately since U is well-constrained. An alternate choice 
is to parameterize the solutions using the dimension¬ 
less Reynolds numbers Ra = aoh/P, Rq = —qnh^//3, 
Rjj = Uqh/P, where ag is the amplitude of Ok. 

For numerical solutions, we adopt the fairly standard 
form 

Ok = ag sin(7rz/ft.). 


We further assume Tc = r, so that 


= 


and take (iKrause fc Raedled 119801 iRuzmaikin et akl 


_ __. 

Brandenburg et al.ll2013|l 


Brandenburg fc SubramaniarJ 


ag = ci^n/h, 


llMl 

( 8 ) 


where C is a constant close to unity; below we set C = 1. 
It is easy to transform to Ra-Rn~Ru parameter space 
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using the relations 


H = 




F = iRr 



but below we work in H-Qt-V space. An important 
quantity is the dynamo number, defined as 

D = RcRn = -q{'ihVL/uf = -q{iHVLTf. ( 9 ) 

If the a term in the dB^/dt equation can be ne¬ 
glected (the a-Q approximation), solutions in the lin- 
ear regime are governe d by the parameters D and Ru 
(|Ruzmaikin et ai.lll988D . 

It is possible to parameterize equations © and © 
in terms of our chosen parameters. To see this, write 
lengths in terms of h (^. z = z/h) and magnetic fields 
in terms of Bq {B = B/Bq), and define = «//? to 
obtain 


dB 1 ~ 


(yzz + Hflrr^) X B 




/ani\ 
dt \ u ) 


2 

3 t 

1 

Hr 


^ siniirl) + ^)b^- xB B 
H u ) 


V • 


3H 

(Viz -k HVlTy)— - ^ 

u aH \ u / 


exp(z 




Thus, if we are interested in averages of B over the 
disk cross section, the solutions are completely deter¬ 
mined if H, T, V, n, q, and R^ are specified. For the 
reader’s convenience, we list in Table [1] the key param¬ 
eters of the models. A typical or ‘canonical’ estimate , 
e.g. for the sol ar neighborhood (|Ruzmaikin et all 119881 
ISliukuro'^l2004n . is given in the rightmost column. 

An asymptotic steady state solution for the sat¬ 
urated regime can be obtained algebraically by ne¬ 
glecting the a term in the (/)-component of Equa¬ 
tion © (the a-II approximation), replacing z-derivatives 
by divis ions by h (with sui t able n umerical coef¬ 
ficients; iSubrani anian fc Mestell Il993t iPhilliosI 120011 : 
iChamandv et al.ll2014a[l . as well as setting time deriva¬ 
tives to zero. Alternatively, the equations can be solved 
algebraically_in the kinematic (linear) regime by writing 
dB/dt = yB, where 7 is the exponential growth rate, to 
obtain 7 and the kinematic magnetic pitch angle pk- Be¬ 
low we sum marize the asymp t otic so lution; details can 
be found in IChamandv et "n] (j2014a[ l. This solution is 
useful as a first line of attack and to gain insight into 
the problem; however, as explained below, it becomes 
inaccurate for some regions of parameter space. 


2 . 1 . Asymptotic Solution 

We now summarize the key expressions from the ana¬ 
lytical no-z solution of iChamandv et ^ (|2014ar) ._J3elow. 
quantities relating to the magnetic held, e.g. p, B^, and 
Bcj,, refer to the steady state (saturated) values, whereas 


the magnetic pitch angle in the kinematic regime is de¬ 
noted by pk- The no-z solution introduces a numeri¬ 
cal coefficient_Ci^_of t he ady ective term in Equation © . 
In IChamandv et HI (j2014a[l . Cu = 1/4 was used, but 
Cu = 1/2 gives a much better ht to the numerically 
determined pitch angles (see Section 0]) , and so we use 
Cu = 1/2 when evaluating expressions. Growing solu¬ 
tions exist if D is larger in magnitude than the critical 
dynamo number, 

D. = -^{'k^ + ^CuRuf = + l2CuHvf . 

32 32 

The magnetic pitch angle p in the steady (saturated) 
state is given by 


1 / tv^ + 12CuHV 

^ Ik ~) ~ 12qH^nT 


( 11 ) 


Note that outhows lead to an increase in |p| because they 
suppress dynamo action, causing \Dc\ to be higher. Note 
also that we have not made explicit use of the dynamical 
quenching nonlinearity © , which may be invoked to ob¬ 
tain an expression for the steady state field strength B, 
which is not needed in the present work. Thus, solutions 
for p are not sensitive to the precise f orm of the dynamo 
nonlinearity (jChamandv et alll2014ajl . 

The growth rate of the dynamo in the kinematic regime 
is given by 


7 = yi 11' (v^ - V^) 

_ ^ [1 + 12CuHV \' 

\ TT J H \288g/ \ ) 

( 12 ) 


where td = /P = 3 H^t is the vertical diffusion time. 
Note that B < 0 for g > 0, that the dynamo is super¬ 
critical for D/Dc > 1, and that larger \D\ corresponds 
to stronger dynamo action. The pitch angle pk in the 
kinematic regime is given by 


tanpk 


ttRq ) 



( 13 ) 


Dividing the first equality of (fOl) by the first equality 
of (fTTI) . and using the definition of the dynamo num¬ 
ber (0 gives 


tanpk / 
tanp \DcJ 


(14) 


so p is similar in both regimes for a mildly critical dy¬ 
namo. 


3. THEORETICAL CONSTRAINTS 

Are there additional constraints on the dynamo equa¬ 
tions, and to what extent do such constraints restrict the 
available parameter space? Eirst, as noted above, the 
dynamo number must have a magnitude that is greater 
than or equal to the critical value, 

\D\ > m, (15) 

otherwise, the magnetic field decays in the linear regime. 
By substituting inequality (fl^ into Equation (fT4l) and 
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Table 1 

List of Key Parameters. Note. Typical estimates, sometimes referred to as “canonical” in the text, refer to the solar neighborhood. 


Description 

Symbol 

Units 

Expression 

Typical Estimate 

Equipartition strength at the midplane 

Bo 

mG 

- 

- 

Magnitude of mean angular velocity 

Q 

km s“^ kpc“^ 

- 

30 

Radial shear parameter 

9 

- 

—9 In r2/91nr 

1 

Mean vertical velocity at the disk surface 

Uo 

km s“^ 


1 

Disk half-thickness 

h 

kpc 

- 

0.5 

Rms velocity of the turbulence 

u 

km s“^ 

- 

10 

Dynamo number 

D 

- 

—^(SD/i/u)^ 

-20 

Correlation length of the turbulence 

1 

kpc 

- 

0.1 

Eddy turnover time 

T 

Myr 

Iju 

10 

Turbulent diffusivity of B 

13 

10^® cm^s-l 

^lu 

1 

Turbulent diffusivity of am 

K 

10^® cm^ 


1 

Dimensionless disk half-thickness 

H 

- 

h/l 

5 

Inverse Rossby number 

- 

- 

Dr 

0.3 

Dimensionless vertical velocity amplitude 

V 

- 

Uolu 

0.1 

Ratio of turbulent diffusivities 


- 

k//3 

1 



Figure 1. Shaded regions denote th e 2 - dime nsi onal projections of the 3-dimensional parameter space that is permitted for growing 
solutions under theoretical constraints n^-iiTsii or 021- (121^ . For ease of presentation, q = 1 is assumed. Dotted lines show how the overall 
upper limits in panels (a) and (b) become modified when K = 1/3 instead of 1. The allowed pa ram eter space under t he m ore restrictive 
conditions with K = 1/3 is shaded darker. Panel (a): Solid lines are from the double inequality (12511 . The lower limit (|23ll becomes more 
restrictive for V > 0. Slices sh own are for V = 0.2 (short dashed), V = 0.6 (long das hed) and V = 1 (dashed-dotted). Panel (b): The 
solid line is from the inequality ll28t . The overall upper limit is equivalent to inequality 02711 with Dr = H. The upper limit becomes more 
restrictive for other values of Dr: Dr = 0.4 (short dashed). Dr = 0.7 (long das hed) , and Dr = 1 (dashed-dotted). Panel (c): the solid 
line is from the inequality ll29t . The upper limit is in general given by inequality 027il . and becomes more restrictive for finite values of H\ 
77 = 5 (short dashed), H = S (long dashed), and 77 = 5/3 (dashed-dotted). Colour figures are available in the electronic version. 


making use of Equation m, we find 


120131 : lOiisholm et al.ll20i^ . and so we require 


tan Ip I < tan |pk| 



(16) 


Relation (fTO tells us that |p| in the saturated regime 
cannot exceed its value in the kinematic regime (lElstnerl 
120051 : l^amandv et al.ll2014all . Thus, although outflows 
help to increasejpl, the latter is limited by |pk|, which is 
independent of Uz in the asymptotic solution. 

Second, the energy density of the helical part of the 
turbulence cannot exceed the total energy density. 


Qfk < Ku, (17) 

where K is a constant of order unity (iMoffattl 119781 
Ruzmaikin et all 119881 iBrandenbure fc SnbramanianI 
200,511 . This' is somet imes referred to as the realizabil¬ 
ity condition (but see iRuzmaikin et al.lll98^ . Below we 
retain K in the equations but for numerical examples we 
adopt A' = 1. Based on equation ((4|), one might argue 
for the stronger condition with K = 1/3, and in some 
figures we show this case for comparison, but we choose 
to be liberal regarding the allowed parameter space given 
the approximate nature of the model. 

Third, we expect the mean vertical veloc- 
ity to correspond t o outflow r ather than inflow 
((Shukurov et al.l 120061 iGent et all l2013al iLaeos et al.l 


Uz > 0. (18) 

We are aware that downward turbulent pumping of B 

would lead to a term with the same form as a down¬ 
ward advective term in Equation Q , and may t hus have 
an im portant effect on the magnetic pitch angle (|Gressell 
[2013). Moreover, some galaxies, especially at high red- 
shift, might accrete, but such complications are left for 
future work. 

If conditions (Ha, (El) and (ITSl) are expressed in terms 
of the dimensionless parameters of interest, we have 

q{mnTf>^{TT^ + i2CuHv)\ (19) 

KH > Hr, (20) 


V > 0. 


( 21 ) 


Below, we also implicitly make use of the facts that H > 
0 and T > 0 and assume q> 0. 

It is useful to rewrite these theoretical constraints as 
constraints on individual parameters. Solving for H from 
relation m, we obtain a lower limit for H valid for all 
rir and V, 


H > 




1/2 


Q.T — 



288q 


( 22 ) 
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On the other hand, obtaining a lower limit for fir 
from inequality CH), and combining this with the up¬ 
per limit (l 20 l) . we obtain 


=; \ 1/2 


1 

li 


12CuV 


<nT< KH. 


(23) 


Given that the left hand side must be less than or equal 
to the right hand side, we deduce a lower limit on H in 
terms of V only: 


1/2 


1 




{CuVf + y(27r3g)l/2 


1/2' 


(24) 

To obtain an overall lower limit on the inverse Rossby 
number Or, we can incorporate condition (EU into ex¬ 
pression (|23|l . so that for all R > 0, we have 


288(/ 


1/2 


— < riT < KH. 
H - 


(25) 


From double 
limit on H, 


inequality (1251) we also obtain an overall 


H > 


tt" 1 

288 gy 


(26) 


where the right hand side evaluates to 1.0 for a flat 
rotation curve (g = 1) with K = 1. This result 
is consistent with the expectation that the turbulent 
scale I should not exceed ^ h for 3-dimensional turbu¬ 
lence. Combining constraint (1^51) with inequality (I16|) 
leads directly to an upper limit to the pitch angle, 
tan|p| < [( 9 /(/)( 2 / 7 r)^]^/^iF^/^ ~ 0.79 for q = 1 and 
iF = 1, or ~ 38°. For q = 0.6, this u pper limit rises 
to^42°. Of the 12 galaxies listed in the iVan Eck et ahl 
(|2015D data, the only galaxy observed to have \p\ larger 
than the above limit is M33, which has p ^ —(36-60)° 
(where the range includes radial variations as well a s 
Icr errors) (|Tabatabaei et al.ll2008l : IVan Eck et al.ll2015l) . 
This galaxy also has rather small angular velocity ft ^ 
25 — 41kms~^kpc~^ and radial shea r q ^ 0.6 — 0.8 
(jSofue et al.lll999l : I Van Eck et al .11201^ . 

Double inequality (1^ governs the allowed H-flr pa¬ 
rameter space, which is a 2-dimensional projection of the 
3-dimensional H-flr-V parameter space. This param¬ 
eter space is represented by the shaded region in Fig¬ 
ure [Ik, where q = 1 and K = 1 have been adopted. A 
dotted line is also included to show how the overall up¬ 
per limit iy = 0) shifts if AT = 1/3 is adopted rather 
than the less conservative A" = 1. The accessible pa¬ 
rameter space under the more restrictive value AT = 1/3 
is denoted by darker shading. If V is specified (and is 
other than 0), then relation (l23)) must be used, and the 
lower limit gets adjusted upwards, as shown by the short 
dashed {V = 0.2), long dashed {V = 0.6), and dashed- 
dotted {V = 1) curves. These represent slices of the 
3-dimensional parameter space. 

We now derive constraints on the velocity ratio V. We 
obtain an upper limit by rearranging expression (HU, and 
then combine this with the lower limit m to obtain 


0 < V < 



Or 


(27) 


Making use of relation (1201) . the upper limit can be gen¬ 
eralized to apply for all r: 


0 < V < 


2gy^^ KH 


12CuH' 


(28) 


This allowed V-H parameter space (for q = 1 and AT = 
1) is shown shaded in Figure [Tla. Expression (l27ll can 
be used to draw slices of this parameter space, that is 
to adjust the upper limit downwards for specific values 
of Or < H. These are represented by the short dashed 
(Or = 0.4), long dashed (Or = 0.7), and dashed-dotted 
(Or = 1) curves. 

Finally, we consider the allowed region of V-Or space. 
For this purpose we require an upper limit on H, which, 
to be liberal, is simply taken to be oo (in practice 
we are more interested in small H anyway). Thus, 
—'k^/{\2CuH) < 0 in the right hand side of expres¬ 
sion (EZl), so that for all H we have. 


0 < V < 


/ 2q\ Or 


(29) 


Thus any upper limit on r translates to an upper limit on 
V. In Figure [!}:, the allowed region of V-Or parameter 
space has been shaded. The solid line corresponds to an 
upper limit on V for H ^ oo, while revised upper limits 
for a sample of finite values of H are also plotted: H = 
5 (short dashed), H = 3 (long dashed), and H = 5/3 
(dashed-dotted). 

The physical constraints (USD, (HZD, and m could be 
recast in terms of p, using Equation (ED, and then ob¬ 
servations of ft, q, and p used to constrain the H~t-V 
parameter space for each data point. However, we leave 
a more detailed comparison of theoretical solutions and 
observational data for a future work. We now turn to 
the main results dealing with theoretical predictions for 
the magnetic pitch angle. 


4. PREDICTIONS FOR THE MAGNETIC PITCH 
ANGLE 

We employ numerical as well as analytic methods and 
compare the values of p from the resulting solutions. 
Closed-form asymptotic analytical solutions can be use¬ 
ful, but it is important to be aware of their limitations, 
and numerical solutions provide an important check. 
Moreover, as mentioned above, the parameter Cjj in the 
asymptotic solution was calibrated using the numerical 
solution. As discussed below, for certain regions of pa¬ 
rameter space the two types of solution differ because 
certain approximations that go into the analytic method 
become inaccurate, but generally, the two types of solu¬ 
tion are mutually consistent so far as predictions of the 
pitch angle are concerned. One minor difference that 
deserves mentioning is that for numerical solutions, con¬ 
straint (I15D is imposed only implicitly, as the field will 
decay if \D\ < |A>c|, while for analytical solutions, it is 
imposed explicitly. In practice, Dc for the numerical so¬ 
lutions will differ slightly from the analytic estimate m 
and must be calculated by iterating until 7 Ri 0. On the 
other hand, constraints (ED/dinD and dUD/dSID are ap¬ 
plied explicitly in the numerical model as in the analytic 
model, that is numerical solutions are not calculated if 
these constraints are not satisfied. 


TT 
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Figure 2. Magnetic pitch angle magnitude, shown as a function 
of the parameters H and Qr, for asymptotic solution (top row) 
and numerical solutions (bottom row), for the case g = 1, and for 
V = 0 (left column), and V = 1 (right column). Parameter space 
for which \p\ > 42° [left side of panels (c) and (d)] is colored black. 
Contours denote 1000-folding times, in galactic rotation periods 
T = 27r/fl, in increments of 2T with dotted curves corresponding 
to 4T and solid curves corresponding to 14T. The thick solid line 
separates the regions for which the kinematic solution is steady 
(below this line) or oscillatory (above this line). Uncolored regions 
designate disallowed parameter space. 
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Figure 3. As Figure[2]but for q = 0.6. 
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Pitch angles are_o_^taiM^ from numerical solutions us¬ 
ing the estimate (pi? )/{B ), where angular brackets rep¬ 
resent the mean value of grid po_u^s between —h and 
h. We weight the pitch angle by B because the inten- 
of polarized synchrotron emission P also scales as 
B . We note that this represents a conservative esti¬ 
mate, as IpI at the midplane is typically larger than_|his 
average value. In fact, to be precise, P (x ncrBj_ds 
where ricr is the number density of relativistic (cosmic 


ray) electron, L is the path length from source to ob¬ 
server, a nd B \ is the co mponent of B in the plane of the 
sky ('e.g. |Shukurc )'\^l2004f) Q If we make the approximation 
(iBeck et al .11 199611 that ricr oc i?^, where B is the strength 
of the tot al (mean-|-random) magnetic field, and assume 
B <x B^r, (iMcKee fc Zwe 1^1199,^ iFletcher fc ShnknrovI 
|200H iBasu fc Rovll2013t iGent et al.ll2013T see also the 
Johns Hopkins Turbulence Database. iLi et ^ 1200811 . 
then the pitch angle should be weighted by the extra 
factor i?gq. We have redone the calculations including 
this extra factor but find that \p\ increases by up to only 
10%, or at most 3° for the parameter space considered. 
In the discussion below, we adopt the more conservative 
flat (z-independent) profile for ricr- 

Figure [2] shows color plots of \p\ as a function of H 
and inverse Rossby number Dr, from analytical (Equa¬ 
tion (HU), top panels (a) and (b)) and numerical solu¬ 
tions (bottom panels (c) and (d)), for the case g = 1. 
Numerical solutions were calculated at steps of 0.01 in 
H and in Dr, with each point corresponding to a sepa¬ 
rate dynamo run with 21-41 grid points (excluding ghost 
zones), as needed. A subset of solutions were checked for 
consistency against higher resolution runs. As expected, 
numerical solutions for p in the steady state, as well as 
for pk in the linear regime, are insensitive to the value 
of which is generally set to 1. The slice V = 0 
is shown in panel (a) or (c) while V = 1 is shown in 
panel (b) or (d). Numerical solutions are generally well- 
approximated by asymptotic solutions, but it can be seen 
that asymptotic solutions tend to underestimate \p\ and 
overestimate \Dc\. Further, the level of agreement be¬ 
tween the two types of solution decreases as V increases. 
From the figure it is evident that theoretical pitch angles 
of magnitude > 20° can be attained, but require values of 
H smaller than canonical, i.e. smaller than ~ 4-5. More¬ 
over, as seen by comparing Figures[2K and[2}D, \p\ is larger 
for given values of PI and r when V is larger. However, 
clearly Dr must be large enough to ensure that the dy¬ 
namo remains supercritical, a constraint that is more eas¬ 
ily satisfied for smaller V. We emphasize that since D is 
generally well constrained by observations, it makes sense 
to discuss the allowed range of r. For example, if Dr = 1, 
then for D = 40kms“^ kpc“^, r = 0.024 Gyr, which is 
somewhat larger than r = 0.01 Gyr, obtained from the 
standard estimates I = O.lkpc and u = 10kms“^. 

We also overplot contours for the 1000-folding time 
Gooo (time for the magnetic field to amplify by the fac¬ 
tor 10^ in the kinematic regime) in units of the galac¬ 
tic rotation period T = 27r/D. Dotted curves corre¬ 
spond to tiooo = 4r and contours are spaced by 2T. 
It can be seen that growth times are over-estimated 
by analytical solutions, especially for V > 0. The 
thick solid line in panel (c) is the dividing line between 
two different regimes. Above this line, kinematic so¬ 
lutions are oscillatory, but become steady upon satu¬ 
ration. This oscillatory behavior is not dependent on 
the effect nor is it sensitive to the for m of aw; the 
fact t h at it was not seen in solutions of [ Brandenbura 
(j 19981 1: [Brandenburg fc Subramani^ (j2005l l is probablv 
due to the symmetry conditions imposed at the midplane 

^ This equation ignores the contribution of the anisotropic ran¬ 
dom component of the magnetic field, but this does not have any 
bearing on the present discussion. 
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in their model, which are different from the ones which 
arise naturally in our solutions. In any case, this regime 
is not particularly relevant for the present work, but its 
investigation in a future study would be interesting. 

The case q = 0.6 is illustrated in Figure|3l Pitch angles 
are somewhat larger in magnitude compared to q = 1. 
This is expected because the magnitude of the fl-effect 
is reduced compared to that of the a-effect. However, 
numerical solutions have noticeably larger \p\ than ana¬ 
lytical solutions for ^ 1 — 2 . As discussed below, this 
is mainly due to the effect, which is neglected in the 
analytic model. 

Figure 0] compares analytical and numerical solutions 
for four sets of parameters, with each set represented 
with a different color. In both panels, curves cut off 
for the value of V above which the dynamo is subcriti- 
cal (where it can be seen that p = pk)- All curves as¬ 
sume H = 40kms~^kpc~^ (close to the median value 
in the iVan Eck et al.l I 2 M 5 I data set). In panel (a), 
solid curves show full steady-state numerical solutions, 
while dashed curves show the corresponding kinematic 
solutions. Occupying the lower portion of the panel 
in blue, we have curves for ‘canonical’ values H = 5, 
T = O.Olskpckm”^ ~ 0.01 Gyr and q = 1. Mov¬ 
ing up to the set of orange curves, we have 77 = 5, 
T = 0.01 skpckm”^ and q = 0.6. The green set of curves 
illustrates H = 5/3, r = 0.03skpckm“^ ~ 0.03 Gyr, 
q = 1. Finally, the upper red set of curves shows 
H = 5/3, r = 0.03skpckm“^, q = 0.6. The simplest 
way to go from the blue to green or orange to red curves 
is to take I to be larger by a factor of 3. Figure HJo 
shows solutions for which the a term was excluded from 
the (/-component of Equation ([3]) (the approxima¬ 
tion), Solid and dashed curves again represent numeri¬ 
cal saturated and kinematic solutions, respectively. We 
also show analytical solutions with Cu = 1/2 (the value 
adopted in calculations, thick dotted), and Cy = 1/4 
(thin dotted). It is apparent that using Cu = 1/2 leads 
to a much more accurate estimate of p. 

Clearly the ra-H approximation is justified for canoni¬ 
cal parameters (blue curves), and the analytical solution 
(with Cu = 1 / 2 ) reproduces well the numerical solution, 
except that it underestimates the maximum V for which 
the dynamo remains supercritical (it overestimates |77c| 
and underestimates 7 ). Going now to the green curves, 
which have smaller H than the blue curves, we see that 
IpI is larger, as expected. The asymptotic solution ap¬ 
proximates quite well the numerical solution that does 
not include the effect, but clearly the (a-H approxi¬ 
mation is inaccurate for large V . Including it leads to 
much larger |p|. (This arises because this term leads to 
a suppression oi near the midplane 0 Finally, mov¬ 
ing to the upper red curves with q = 0 . 6 , it is evident 
that IpI is larger than for g = 1, as expected. More¬ 
over, the c? effect plays an even bigger role since the 
ratio _of the shear term to the a term in the equation 
for dB^/dt is proportional to q. We note that if the 
dynamo is near-critical, the growth rate is small, and, 
depending on the seed field assumed, the magnetic field 

^ In principle this effect could be accommodated in the analytic 
model by replacing q with q — g, where g is a positive function of 
some of the parameters. However, the nonlinear dependence of p 
on V complicates such an approach. 


may not have time to saturate for a given integration 
time (say 10 Gyr). However, this does not mean that 
large magnitude pitch angles cannot be achieved in such 
cases, because |p| approaches monotonically or almost 
monotonically the saturated value from its value in the 
kinematic regime |pk|, which is, in any case, the upper 
limit of |p|. 

The spatial profile of the magnetic field for each of the 
parameter sets of Figure |4] is plotted in Figure [5l us¬ 
ing the same color for each paramet^set asjn Figure SI 
Gurves for the steady state values of B^ and 77^ are plot¬ 
ted on the left, while p is plotted on the right, for H = 0 
(solid), V = 0.1 (short dashed), V = 0.2 (dash-triple- 
dottecl), V = 0.4 (long dashed), and V = 0.8 (dashed- 
dotted), for cases for which such an outflow does not 
lead to a subcritical dynamo. All quantities_plotted are 
even about the midplane. In om solutions, 77^ tends to 
mainly negative values, while 77^ tends to mainly posi¬ 
tive values, but the equations are invariant under a sign 
reversal of B, and the signs ar^determined by the ar¬ 
bitrary seed field chosen. Both 77^ and 77^ are normal¬ 
ized with^espect to the midplane mean field strength 
77(0^ =J^r( 0 ) + The average pitch angle 

(p77 )/(77 ) is also shown as a horizontal line of the ap¬ 
propriate linestyle. As mentioned above, the largest val¬ 
ues of IpI arise for small H and large Hr and H, when 
there is a significant reduction of B^ near the midplane, 
a consequence of the o? effect. 

We see then that large magnitude pitch angles can in¬ 
deed be achieved for small values of 77 and large values 
of Hr and V relative to canonical. This cannot be fully 
appreciated by using the asymptotic solution because it 
does not account for the effect, and so underesti¬ 
mates IpI, especially for the parameter values that are 
favourable for large pitch angles. Moreirver, p ast works 
(|Chamandv et al.ll2014at IVan Eck et al.l 1201^ tend to 
underestimate |p| for non-zero outflows because they use 
Cu = 1/4, whereas here we use Cu = 1/2, which is a 
much better calibration to numerical solutions for mak¬ 
ing estimates of p (though perhaps worse for making es¬ 
timates of 7 ). Perhaps the simplest way to achieve large 
IpI is to make I and Uq larger than standard estimates. 
Is such a prescription plausible in light of independent 
evidence? 

5. WHAT PARAMETER VALUES ARE REALISTIG? 

The magnetic field is believ ed to roughly tr ace the 
warm ionized gas component (iBeck et al.l Il996ll . Out¬ 
flow speeds of the hot gas tond to be as large as a 
fewxl0^kms“^. However, Uz represents an average 
over all interstellar medium (ISM) phases inhabited by 
the mean magnetic field, as well as over the disk sur¬ 
face, which includes regions in between hot sup^bub- 
bles (chimne ys), leading to standard estimates 7/„ ^ 
0.2-2 kms~^ (lShukurov et al.ll2006l : IVan Eck et al.ll2015l : 

iRodrigues et al.ll2015[l . Such estimates involve equating 
the mass flux for the averaged ISM to that of the out¬ 
flowing hot gas. However, simulations and observations 
seem to indicate that warm gas can become entrained in 
the hot outflows, leading to velocities of the warm gas of 
several tens of kms~ ^ or more (e.g. iGent e^l.l[2013al : 
iGhisholm et al. I [20141) . Therefore, values of Uz compa¬ 
rable with u ^ 10-20 km s“^ could possibly be more re¬ 
alistic than past estimates, at least for some galaxies. 





















Figure 4. Panel (a): Full numerical saturated solution 

(solid), as well as full numerical kinematic solution (dashed), for 
four different sets of parameters, distinguished by color as indi¬ 
cated on the plot. All curves have Q = 40kms“^kpc“^ and 
= 1 (though solutions are insensitive to Rk)- Lines are drawn 
from V = 0 up until the point where the dynamo becomes critical. 
Panel (b): Numerical saturated solution with approximation 
(solid), along with the corresponding kinematic solution (dashed). 
Analytical steady state solutions are shown for Cu = 1/2 (thick 
dotted) and Cjj = 1/4 (thin dotted). For both values of C[/, |p| 
rises to the same value |pk| before the solution goes subcritical. 


It is worth emphasizing, however, that an extra term of 
the form 7 x B in Equation ([2]) for with Iz < 0 , 
could b ecome impor tant when the gaseous halo is in¬ 
cluded (|Gresselll2010ll . The resulting turbulent pumping 
of the magnetic field towards the midplane would lead 
to an effective value of Uz in Equation dS]) that is lower 
than the actual mean vertical velocity, and could, in prin¬ 
ciple, even be negative. This pumping term would not, 
however, play any role in Equation ([S|), so the advective 
helicity flux would not be directly affected. 

The correlation length of turbulence I has been rela¬ 
tively challenging to determine observationally. Many 
studies have used Faraday rotation measurements to 
probe the scale of magnetic field fluctuations. How- 
ever, these scales may be differen t from the scale I 
(|Brandenburg fc Subramanianl 1200^ . More relevant to 
the present study, perhaps, are studies that set out to 
determine the ve locity correlation l ength directly. For 
the Milky Way, iChepurnov et al.l (j201(lH find a scale 
of 0.140 ± 0.080 kpc for the warm HI. Observations of 
nearby dwarf galaxies yield values for I of a few kpc 
(lElmegreen et all 120011 : Stanimirovic fc LazarianI 120011 : 
iDib fc Burkerd |2005[) . Dutta et al.l (|2013li obtain HI 
power spectra for a sample of 18 spiral galaxies and find 



Figure 5. Magnetic field profiles (left) and magnetic pitch angles 
p (right) for each of the parameter sets depicted in Figure |4l with 
corresponding color. Panel (a)-. H = 5/3, r = 0.03skpckm“^, 
q = 0.6; Panel (b): H = 5/3, t = 0.03skpckm“^, q = 1; 
Panel (c)\ H = 5, t = 0.01 skpckm~^, q = 0.6; Panel (d): 
H = 5, T = 0.01 skpc km“^, q = 1. Linestyles represent different 
outflow speeds: V = 0 (solid), V = 0.1 (short dashed), V = 0.2 
(dash-triple-dotted), V = 0.4 (long_d|shed), and V = 0.8 (dashed- 
dotted). Horizontal lines show the B -weighted spatial mean of p. 


power law fits extending up to a scale comparable to the 
radial scale length of the optical disk (several kpc), with 
power laws flatter than those found at smaller scales. 
They attribute this flattening to a transition from 3- 
dimensional to 2 -dimensional turbulence for I ~ h, where 
h is the disk scale height. 

Many different drivers of tur bulence in the ISM have 
been proposed in the literature (|Falceta-Goncalves et al.l 
[Ml and references therein), and more than one such 
driver may, in fact, be important. Hydrodynamic sim¬ 
ulations of the local ISM with supernova-d riven tur¬ 
bulence by Ide Avillez fc Breitschwerdtl (I2007D pre dict a 
value I ~ 0.075 kpc. whereas iGent et al.1 ( 20i3a[) pre¬ 
dicts I ^ 0.1 kpc at the midplane, rising to ^ 0.3 kpc at 
1^1 = 1 kpc. These determinations of I were for the aver¬ 
aged ISM in the solar neighborhood, rather than for any 
particular phase. Simple analytic estimates based on the 
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size attained by a supernova remnant by the time it mixes 
with th e ambient IS M give I ^ 0.1-0.2kpc for the warm 
gas fe.g. lLacl^l2013D . However, if the turbulence is driven 
by expanding supernova-driven superbubbles, then the 
turbulent scale may be comparable to the largest sizes 
attained by such bubbles before mixing with the ambi¬ 
ent ISM. Such bubbles may become as l arge a s a few kpc 
in some galaxies (e.g. iBoomsma et al.l 1200811 . exceeding 
the scale height of the thin component of the disk, but 
still comparable to t hat of the diffuse ionized medium 
(|Gaensler et all I2008H . The bubble cross-sectional area 
at the mixing stage is expected to increase with distance 
z from t he midplane, a s well as with the galactocentric 
radius r (lFerrierelll998fl . It is not clear whether the box 
size of the ISM simulations mentioned would have been 
large enough to ca pture turbulence on suc h large scales. 
In a recent work, iMoraehan et al.l (|2015ll find in their 
simulations that the spherical outflows that drive the 
turbulence transfer power to scales larger than the initial 
injection scale as they expand. Based on all of this evi¬ 
dence, it is not possible to exclude any value ioi H = h/l 
in the range I-IO. 

Note that raising I has no effect on the dynamo num¬ 
ber D (Equation dH])), though it does have the effect of 
decreasing the magnitude of the critical dynamo number 
Dc (Equation (ITUll b Thus, the ratio D/Dc, which gov¬ 
erns dynamo growth via Equation (1121) . increases from its 
canonical value, if other parameters are kept constant. In 
fact, it can be seen from Equation (TT^ that the exponen¬ 
tial growth rate of the mean magnetic field in the kine¬ 
matic regime 7 increases for two reasons. First, \Dc\ is re¬ 
duced, and second, the turbulent diffusion time td o^ 1/1 
is reduced. Thus, a larger I implies a more supercritical 
dynamo which is able to grow the mean magnetic field 
more rapidly. This is also evident in Figures [2] and [3l 
where smaller H and larger r correspond to smaller ex¬ 
ponential growth times in the kinematic regime. On the 
other hand, a model with a larger I can accommodate 
a larger, perhaps more realistic value of E, for a given 
growth rate. 

The tentative conclusion is that the parameter values 
required for larger magnitude pitch angles are not in¬ 
consistent with observational and theoretical constraints. 
However, these parameter values (most notably H ^ 
1 — 3) sit in a small region of the allowed parameter space, 
near a ‘corner’ (left-hand sides of the panels of Figures [5] 
and El) that is formed by rather robust theoretical con¬ 
straints (‘slanted walls’), namely the requirement for dy¬ 
namo growth and the realizability condition. Therefore, 
although IpI > 20°, and even \p\ > 40° as in the galaxy 
M33, is in principle attainable, the parameter values that 
allow large pitch angle are at the extreme range of the 
allowed parameter space. For the largest possible val¬ 
ues of IpI, the dynamo must be close to critical and the 
turbulent velocity field must be nearly maximally helical 
(even then if if = 1/3 instead of if = 1, \p\ > 40° would 
no longer be attainable). Of course, it is possible that 
physical mechanisms exist to push parameters into this 
region of parameter space. It is also possible that the 
theory is still missing some important physical ingredi¬ 
ents that would lead to new effects that would in turn 
tend to increase the value of \p\. Where should we look? 

6 . POSSIBLE EXTENSIONS OF THE MODEL 


There are many possibilities for extending our basic 
model. The most obvious is to include the radial de¬ 
pendence. This would introduce new parameters (radial 
scale lengths) that would complicate the model. Past 
studies that included the radial coordinate have shown p 
to be of roug hly the same magnitude as predicted by lo¬ 
cal solutions (iChamandv et al1l2013ll . To verify this, we 
have carried out a preliminary study using a mean-field 
simulation in 1 —‘no- 2 :’, with a Brandt rotation curve (see 
IChamandv et al.ll20I3l. for details of the model). Outside 
a small central region (where the thin disk approxima¬ 
tion may not be valid anyway) we generally find that 
p{r) in the saturated regime is very close to the local an¬ 
alytical prediction, until r > Tc, where Tc is defined by 
D{rc)/Dc{rc) = I. For values of r up to ^ I kpc larger 
than Tc, the magnetic field in the saturated state can 
still be significant, and interestingly, we find that here 
the pitch angle saturates near its kinematic value. While 
such issues are beyond the scope of this paper, they de¬ 
serve further study. Specifically, detailed comparisons of 
solutions obtained for local models in z and global mod¬ 
els in 1 —z with and without the thin disk approximation 
would be useful. It would also be interesting to com¬ 
pare theoretical predictions for the variation of p with 
galactocentric radius to present and future observational 
data. 

Another possible extension of the model is to incorpo¬ 
rate a gaseous halo (corona) and turbulent pumping (tur¬ 
bulent diamagnetism). One effect of this would be the 
ability of the d ynamo to withstand larger outflow speeds 
(iGressell[20I0H . The effect of a gaseous halo, including 
turbulent diffusivity /3 that varies with height, has not 
been studied in the context of dynamical quenching the- 
ory. 

In certain ways, our model may be oversimplified. 
For example, the turbulent correlation time may 
not, in fact, be equal to the eddy turnover t ime r for 
supernova-driven turbulence (|Shukurovl l2004f) . More¬ 
over, the parameters H, r, and V may vary signifi¬ 
cantly with position from the midplane z. It also de¬ 
serves to be mentioned that turbulent transport coef¬ 
ficients such as a (or in general aij) are known to 
experience rotational quenching and to become more 
anisotropic as Hr increases beyond I (iRuzmaikin et aH 
119881 [Brandenburg &: Subramanianll2005L and references 
therein). As expressions in our model are, in any case, 
imprecise to within factors of order unity, and as we have 
restricted our analysis to values of Dr < 2 , the neglect of 
these effects is not an important limitation of the present 
study. However, possible modifications to the model for 
large Hr are worth exploring in a future study. 

Introducing non-axisymmetry into the dy¬ 
namo equations may result in a number of ef- 
fects, including new contributions to the £ term 
(iBrandenburg fc SubramanianI 1200511 . In particular, 
spiral shocks wou ld tend to align mag netic field along 
the spiral arms ([Van Eck et al.l I20I5I1 . and generally 
streaming flows associated with spiral density waves 
might lead to important effects on B. It should be 
noted, however, that at least for the most extreme case 
of M33, alignment of magnetic field with spiral arms 
would not help because the pitch angle of the spiral 
arms is smaller in magnitude than that inferred for the 
mean magnetic field of that galaxy ([Van Eck et al.ll2015L 
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and references therein). 

One still highly contentious area of mean-field dy¬ 
namo theory is the nature of the magnetic helicity 
flux and the corresponding flux density term in the 
dynamical quenching Equation m for a^n- As we 
have shown, the value of p is independent of the 
nonli nearity used for a cer tain class of nonlineari¬ 
ties (|Chamandv et al.l l2014all . However, it is pos- 
sible that our model ex cludes important flux terms 
( [Vishnia£_fc_Cl^ [20^ [Subramaniai^_^^randenburd 
20001: iV ishniad l2012t IVishniac fc Shapovalov! 12014 
Ebrahimi fc Bhattacharied I2014fl . that may drastically 
affect the mean magnetic field and its pitch angle. 

In fact, we are aware of a few alternate models that 
potentially lead to large pitch angles. A mean-field dy¬ 
namo for which the a term is dominated by the Vishniac- 
Cho (V C) flux leads to p ^ — tan“^[7r/(2iJ)] ^ —17° for 
H = 5 (|Sur et al.l 1200711 . However, it is not clear how 
such a dynamo w ould saturate. The dynamo model of 
iMoss et (|1999fl based on the buoyancy-driven Parker 
instability results in estimates of the pitch angle that can 
be very large (their Equation (19)). This model is worthy 
of more exploration in the future. Notably, both the VC 
flux and Parker instability mechanisms exhibit threshold 
behavior, in that they require a mean magnetic field of 
a certain strength to be present initially. 

Another promising avenue for producing large pitch 
angles is the magnetorotational instability (MRI) 
(|Kitchatinov fc Riidigeii [200l: lElstner et al.ll2009ll . The 
unstable MRI modes have p = —45° in the ideal MHD 
limit, but including finite viscosity and diffusivity (e.g. 
due to turbulence) might result in considerably smaller 
IpI (|Pessah fc ChanI [200811 . In any case, the influence of 
the MRI on the galactic magnetic field deserves further 
study. 

7. SUMMARY AND CONCLUSIONS 

Magnitudes of magnetic pitch angles of the large-scale 
magnetic fields of spiral galaxies are generally underes¬ 
timated by nonline ar dynamo theory. F or the 31 data 
points compiled bv IVan Eck et al.l (I20I5I1 . p has a mean 
value of —25° and minimum and maximum values of 
—8° and —48°, excluding uncertainties, whereas stan¬ 
dard nonlinear dynamo theory with canonical parame¬ 
ters generally predicts \p\ < 15°. 

We have taken a fresh look at the predictions of p 
from standard mean-field galactic dynamo theory with¬ 
out making any a priori judgements about parameter 
values. Our basic model is a thin disk (slab) dynamo sub¬ 
ject to magnetic helicity balance, that grows the mean 
magnetic field to a steady saturated state. We have 
shown that in this model solutions for p in the saturated 
state (or in the kinematic regime for that matter) can 
be usefully parameterized in terms of the angular veloc¬ 
ity H, shear parameter q = —cHnU/cHnr, scale height 
of the disk in turbulent correlations lengths H = h/l, 
turnover time of energy-carrying eddies r = l/u, and 
mean outflow velocity in units of the rms turbulent ve¬ 
locity of energy-carrying eddies V = Uq/u. Of these, H 
and q are well-constrained observationally, leaving a 3- 
dimensional parameter space, which is dimensionless if 
the inverse Rossby number Hr is used in place of r. The 
range of allowed parameter space is further restricted by 
theoretical considerations, namely the requirements that 


the mean magnetic field grows rather than decays and 
that the kinetic energy in the helical part of the turbu¬ 
lence is less than or equal to the total turbulent kinetic 
energy. 

This still leaves the standard parameter space as well 
as a rather small non-standard region of parameter space 
that leads to \p\ > 20°, and a much smaller region that 
gives \p\ > 40°. To be specific, such large \p\ requires 
smaller than canonical H and larger than canonical r, 
which can be achieved simultaneously and rather natu¬ 
rally by adopting a turbulent correlation scale I that is 
a few times its canonical value of O.Ikpc, if the stan¬ 
dard value h = 0.5 kpc is assumed. Larger than canon¬ 
ical values of V are also required for the largest |p|; on 
the other hand, dynamos with larger I have larger kine¬ 
matic growth rates and are better able to withstand a 
large V without becoming subcritical. We have argued, 
based on independent evidence from the literature, that 
such changes to the parameter values are plausible. How¬ 
ever, the need to adopt parameters near one extreme of 
the allowed parameter space suggests that physical el¬ 
ements may be missing from this basic mean-field dy¬ 
namo model, and we have offered several suggestions for 
extending the model to include potentially important ef¬ 
fects. 

Our results were obtained using asymptotic analytical 
as well as numerical solutions, and are almost indepen¬ 
dent of the precise form of the dynamo nonlinearity as¬ 
sumed. Good agreement was obtained between the two 
types of solution for most of the parameter space. The 
analytical solution becomes inaccurate for parameters for 
which the effect becomes important. Interestingly, 
these are the parameters that lead to the largest values 
of Ip I, rendering numerical solutions indispensable. 

A more detailed comparison between theory and obser¬ 
vation is warranted. A likelihood an a lysis o f the param¬ 
eter space, given the IVan Eck et al.l (|20I5ll data, would 
be interesting both in its own right and as a demonstra¬ 
tion of what could be done with better data sets and 
perhaps better models. A related goal would be model 
comparison between the basic model presented and mod¬ 
els that incorporate new physical effects, but as a result 
may require more free parameters. 
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